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Abstract 



We study Low-Density Parity-Check (LDPC) codes with iterative decoding on block-fading (BF) 
0^ ! Relay Channels. We consider two users that employ coded cooperation, a variant of decode-and-forward 

with a smaller outage probability than the latter. An outage probability analysis for discrete constellations 
shows that full-diversity can be achieved only when the coding rate does not exceed a maximum value 
Q ' that depends on the level of cooperation. We derive a new code structure by extending the previously 

published full-diversity root-LDPC code, designed for the BF point-to-point channel, to exhibit a rate- 
^ . compatibility property which is necessary for coded cooperation. Word error rate (WER) performance 

■ is determined for infinite length codes through density evolution (DE) as well as for finite length codes. 

' We show that our code construction exhibits near-outage Umit performance for all block lengths and for 



a range of coding rates up to 0.5, which is the highest possible coding rate for two cooperating users. 
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I. Introduction 



When communicating over fading channels, Word Error Rate (WER) performances as well 
as power savings are dramatically improved through transmit diversity, i.e., transmitting signals 
carrying the same information over different paths in time, frequency or space. Spatial diversity 
obtained by the transmitter generally requires a Multiple Input (multiple antennas) system. 
However, many wireless devices have constraints in size or hardware complexity so that they 
can be equipped with only one antenna. Recently, a new network protocol called Cooperative 
Communication [11], [22], [29], [36], [37] yields transmit diversity using single-antenna devices 
in a multi-user environment by taking advantage of the broadcast nature of wireless transmission. 
The most elementary example of a cooperative network is the relay channel, introduced by 
van der Meulen [28]. In a relay channel, a relay helps the source in transmitting its data 
to a destination by relaying the messages sent by the source so that the received energy at 
the destination is increased. This relay channel can be generalized to a cooperative Multiple 
Access Channel (MAC) [22], where two users transmitting data to a single receiver cooperate by 
alternately being the relay for the other user, as indicated in Fig.H] Further generalization to more 
users is possible, but this will not be discussed here for simplicity. Cooperative communications 
with practical algorithms in a wireless context was first shown to have a higher throughput than 
direct transmission by Sendonaris et al. [36], [37]. 



Fig. 1. A Cooperative Multiple Access Channel (MAC). Arrows between two nodes illustrate that both nodes communicate 
between each other. 



A challenging channel model is the BF [4] frequency non-selective Single-Input Single-Output 
(SISO) channel. When the fading gain is constant over a codeword and no cooperation is used. 



User 1 




Destination 



User 2 
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the resulting word error rate curve (displaying the logarithm of the error rate versus the average 
signal-to-noise ratio (SNR) in dB) has the same high-SNR slope as for uncoded transmission: 
the corresponding diversity ordeu equals one. In this case, a lot of transmit energy is needed 
to achieve sufficiently low error-rates. Automatic Repeat-reQuest (ARQ) is one of the means to 
increase the diversity, at the expense of increased packet delay. The potential diversity increase 
brought by cooperative techniques allows to save much transmit energy at a given error rate 
and to reduce the delay of the received packet due to less retransmission requests. This channel 
model is adopted here because of the large promised gains due to cooperation. BF channels are 
a realistic model for a number of channels affected by slowly varying fading and flat fading is 
assumed in order to isolate the effect of cooperative diversity. 

The specific task of the relay is determined by the strategy or protocol. In the case of Amplify 
and Forward (AF), the relay amplifies the received signal and noise; for Decode and Forward 
(DF), the relay first decodes and then re-encodes the message before sending it to the destination. 
It is not straightforward to state which strategy is optimal, because this depends on the channel 
conditions, such as the ergodicity of the temporal correlation of the fading gain and the distances 
between the source, the relay and the destination [21]. A variant of DF is coded cooperation, 
where the relay decodes the message received from the source, and then transmits additional 
parity bits of the message, resulting in a more spectral efficient strategy [18], compared to a 
traditional DF protocol. Instead of SNR accumulation (logarithmic rise of mutual information 
with received power from the relay) at the destination, we get information accumulation (linear 
rise of mutual information with received power from the relay) [42]. It has been shown in [19] 
that the outage probability [4], [30] of coded cooperation for half-duplex BF channels is smaller 
than for repetition-based protocols. Moreover, the concept of coded cooperation can be used 
in more complex strategies, such as Amplify-Decode-Forward [2], where the relay can choose 
between DF and AF. So finally, replacing the decode-and-repeat part in any protocol by this 
more intelligent "information adding" strategy improves the outage probability performance. As 
a consequence, constructing a near-outage channel code for a coded cooperation scenario results 

'Here, diversity order is defined as the ratio of the high-SNR slopes of the error rate curves of the considered system and 
of the uncoded system, respectively. Alternatively, diversity order can be defined as the slope of the error-rate curve of the 
considered system. The diversity depends on the fading gain distribution in the latter definition, but not in the former definition. 
Both definitions are equal in the case of Rayleigh fading. 
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in a competitive error-correcting code in terms of error-rate performance vs. SNR for a given 
rate R. 

In a first part, this paper analyzes the outage probability for binary phase shift keying (BPSK) 
modulations and derives a coding rate limitation that is necessary for the protocol to have 
diversity two, valid for all discrete alphabets. But not every code satisfying this information 
theoretical rate limitation achieves a double diversity order. In contrast to DF, the diversity order 
is not guaranteed by the protocol itself. Deriving a code structure exhibiting full-diversity will 
be treated in the second part of the paper. 

Up till now, coded cooperation has mainly been implemented using rate-compatible convo- 
lutional codes [18]. The main drawback of these codes is that the WER increases with the 
logarithm of block length to the power d where d is the diversity order [7], [8]. The WER 
of practical near-outage codes should be independent of the block length in order to approach 
the outage probability limit [12], [13]. The solution is to use capacity-achieving codes, such as 
LDPC codes [32] or Turbo Codes [3]. In this paper, we concentrate on the family of LDPC 
codes. LDPC codes designed for the special case of a cooperative channel have been reported 
by Razaghi et al. [31] for the Gaussian channel. For the block- fading channel however, there is 
still a lack of a near-outage LDPC code. Hu et al. [16] obtained interesting results by designing 
random LDPC codes for cooperation over fading channels, but they did not tune the parity- 
check matrix in order to achieve full-diversity, as shown by Boutros et al. [6] and as will be 
explained in the rest of the paper. The aim of coded cooperation is to send a codeword over two 
independent fading paths and the relay must be able to decode after receiving the first part of 
the codeword. An error-correcting code must therefore exhibit two properties: full-diversity and 
rate-compatibility. This paper derives a new code structure satisfying both properties. We also 
determine density evolution equations to study the performance of an infinite length LDPC code. 
These results are representative for finite length codes [32] as illustrated in the numerical results. 

Channel-State Information (CSI) is assumed at the decoder. We consider half-duplex devices, 
assuming that simultaneously receiving and transmitting data in the same frequency-band is too 
complicated due to the limited isolation of directional couplers. In addition, we also restrict the 
protocol to be orthogonal since we transmit at low rates (we use Binary Phase-Shift Keying 
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(BPSK)). The proposed code construction can nevertheless be used in more complex non- 
orthogonal protocols, where one can achieve more coding gain in high-rate scenarios [1]. 

II. System model and notation 

As mentioned in the introduction, the devices are half-duplex and users transmit in non- 
overlapping time slots. The transmission of a codeword is organized in two frames which 
constitute one block. We denote the transmission of user u, u = 1,2, in frame m, m = 1,2, 
by Xu^m. The pairs {Cu,i,Cu,2) denote the codewords of user u. In the first frame of a block, 
each user broadcasts the first part of its encoded data to the other user and to the destination. 
In the second frame, users either cooperate or send additional parity bits related to their own 
information message, depending on whether they are able to decode the transmissions in the first 
frame. The decoding failure is detected by the relaying user via a Cyclic Redundancy Check 
(CRC) code or any other intelligent detection scheme. There are 4 cases to be distinguished, 
as summarized in Fig. [21 in case 1, both users have successfully decoded the information from 
the other user; in case 2, none of the users has been able to decode the information from the 
other user; in case 3 (case 4), only user 2 (user 1) has successfully decoded the information 
from the other user. Two methods have been proposed in [17] allowing the destination to detect 
which of these 4 cases has occurred:(i) The first method is to decode according to all possible 
assumptions until CRCs indicate successful decoding, (ii) The second method is to let each user 
transmit one additional bit to indicate its state to the receiver. 

A codeword will consequently be split over 2 frames. We consider codewords to have a total 
length equal to N binary digits, where = A'^i + and A^i and N2 denote the length of the 
first and second part of the codeword. We define the level of cooperation, /3, as the ratio N2/N. 
We denote the transmitter of a frame, which can be user 1 or user 2, by s and the receiver of 
a frame, which can be user 1, user 2 or the destination, by r. Transmitted symbols of user 1 
will be denoted xi[i] where i is the symbol time index, i E {1, . . . , A^}. Similarly, transmitted 
symbols of user 2 are denoted X2[i]. The transmitted symbols are chosen from a BPSK alphabet, 
Xs[i] E {1,-1}. Received symbols will be denoted ysr[i] for received symbols from transmitter 
s to receiver r. The received symbol is given by 

ysr[i\ = asrXs\i] + Zr[i\, (1) 
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Frame 1 


Frame 2 




Frame 1 


Frame 2 


User 1 


-^1.1 — C*!,! 








User 1 


^1.1 — C^Ll 




X\.2 — C12 




User 2 




^2.1 — C2.I 




X2.2 = Cl.2 


User 2 




-^^2,1 — C2.\ 




X2,2 — C2:2 



(a) Case 1. Both interuser transmissions are successfully 
decoded. Each user cooperates in the second frame. 



User 1 



User 2 



Frame 1 



-'^i.i — Ci.i 



^2.1 — ^2.1 



Frame 2 



X1.2 = c\. 



X2.2 = Ci, 



(c) Case 3. User2-to-Userl communication failed. In the 
second frame, user 1 sends its own parity bits and user 2 
cooperates witlr user 1. 



(b) Case 2. Both interuser communications failed. Each 
user sends its own parity bits in the second frame. 



User 1 



User 2 



Frame 1 



Frame 2 



^1,1 ~ Ci i 



Xl,2 = C2. 



X2.I = C2. 



X2:2 — C2: 



(d) Case 4. Userl-to-User2 communicatino failed. In the 
second frame, user 2 sends its own parity bits and user 1 
cooperates with user 2. 



Fig. 2. The 4 cases encountered in coded cooperation are listed above. 



where Zr\i\ ~ A/'(0, cr^) are independent noise samples and G is the Rayleigh distributed 
fading gain between sender s and receiver r. The fading coefficient asr is assumed to be constant 
during 2 frames and is independent from block to block due to interleaving. 

We focus on binary LDPC codes C[N, K]2 with block length A^, dimension K, and coding rate 
Rc = K/N . The code C is defined by dM {N — K) N parity-check matrix H, or equivalently, 
by the corresponding Tanner graph [5], [33]. It is assumed that H has full rank (A^ — K), and 
consequently i?c = 1 — ^'n^ ■ Regular LDPC ensembles are characterized by the pair (db,dc), 
where db is the maximum bitnode degree and dc is the maximum checknode degree. L^regularity 
is introduced through the standard polynomials \{x) and p(x) [34]: 

db dc 
i=2 1=2 

where Aj (resp. pi) is the fraction of all edges in the Tanner graph, connected to a bitnode (resp. 
checknode) of degree i. Therefore, A(x) and p{x) are sometimes denoted as left and right degree 
distributions from an edge perspective. In section |V] the polynomials A(x) and p{x), which are 
the left and right distributions from a node perspective, will also be adopted: 

dh dc 
i=2 1=2 
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where Aj (resp. pi) is the fraction of all bit nodes (resp. check nodes) in the Tanner graph of 
degree i, hence = j2 'x./j and likewise with pi. 

In this paper, not all bit nodes and check nodes in the Tanner graph will be treated equally. 
To elucidate the different classes of bit nodes and check nodes, a compact representation of the 
Tanner graph, adopted from [9] and also known as protograph representation [38], will be used. 
In this compact Tanner graph, bit nodes and check nodes of the same class are merged into one 
node. 



The BF channel has a Shannon capacity that is essentially zero since the fading gain makes 
the mutual information a random variable which does not allow us to make the error proba- 
bility arbitrarily small under a certain spectral efficiency. This word error probability is called 
information outage probability in the limit of large block length, denoted by 

Po„i = P(/(a,7)<i?), 

where /(a, 7) is the instantaneous mutual information as a function of a certain fading gain a 
and average SNR 1^ 1 = = where Eg is the symbol energy. This definition remains valid 
for a channel model as described in O, but then a is the set of fading gains over a codeword 
and 7 is the set of average received SNRs. The rate R is the spectral efficiency of a user, only 
taking into account its timeslots, hence not the average spectral efficienc>0. It is our aim in this 
paper to approach the outage probability limit for a range of values of the spectral efficiency R. 
Since we use BPSK signaling, the spectral efficiency R is identical to Re- 

Definition 1 The diversity order attained by a code C is defined as 

a = — iim , 

7— >oo log 7 

where is the word error rate after decoding. 

Definition 2 An error-correcting code is said to have full- diversity if d = N^, where iV„ is the 
number of cooperating users. 

^This is, in our opinion, necessary for a fair comparison between multiple user networks with a different number of users. 
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Notice that the above definition assumes Rayleigh distributed single antenna channels. The 
diversity order of the outage probability limit is the same as the order attained by a full-diversity 
channel code [12]. According to the blockwise Singleton bound [20], [27], the coding rate for 
a ra-order full-diversity code is upper bounded by Rcmax = ^/n. Hence, in a 2-user scenario we 
get Rc < 0.5. 



The outage probability analysis of coded cooperation with a Gaussian alphabet has been 
made in [19]. Here, the analysis considers BPSK signaling, leading to an important conclusion 
in Corollary \T\ at the end of this section. The stated corollary is also valid for larger discrete 
alphabets. 

The average mutual information of a SISO channel with received signal y, conditioned on the 
channel realization a, is determined by the following well-known formula [40]: 



where 'Eyia is the mathematical expectation over Y given a. The outage event of a point-to-point 
link is defined by the mutual information of that link being less than its transmission rate. The 
outage event Eo of the relay channel is determined by a specific region in the multidimensional 
space of instantaneous signal-to-noise ratios. Next, we give the exact definition of Eo for coded 
cooperation with BPSK modulation. We shorten the notation /(X^; Yj\aij) to lij. 

Proposition 1 In coded cooperation for a two-user MAC with BPSK signaling, the outage event 
Eo related to user 1 is expressed as follows: 



III. Outage Probability Analysis 




(2) 



E, 



o 



a 




p))n{hdicase 1)<R)] 



U [(/i2 < R/{1 



P)) n {hi < R/{1 



P))n{hd {case 2)<R)] 



U[(Ji2 >/?/(! 



(3)) n {hi < R/{1 



p)) n {hd {case 3) < R)] 



U[(/i2<i?/(l 



/?)) n (/21 > R/{i 



P))n{hd{case4)<R/{l-m }, (3) 
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where 



I 



b) 



1 - Ey|ai2 <^ log2 1 + exp 



12 



12 



^ 1 
21 — J- 



Ey|a2i < logs 1 + exp 



-2l/21«21 



0". 



21 



(4) 
(5) 



and where lid is case-dependent. For each of the cases considered in Fig. |2] the mutual 
information Iid can be calculated as follows: 



Case 1: 



Ey,,^, aog2 1 + exp 



-2yidaid 



a 



Id 



Case 2: 



-'2yida 



a 



Id 



Jirf = 1 - Ey|„^^ <^ logs 1 + exp 



Case 3: 



/id = 1 - (1 - P) ^Y\a,, < log2 ( 1 + exp 



Id 



Id 



<^ logs 1 + exp 

, (ttldZ/ld + «2dl/2d) 



-2{y')iald + aldr^' 



V"?d + 



a 



2d 



Case 4: 



- (3 Ey|„2rf <^ logs 1 + exp 



-2y2da2d 



2d 



(6) 



(V) 



(8) 



/id = 1 - Ey|„^^ <^ logs 1 + exp 



-2?/ida 



Id 



Id 



(9) 



Proof: 

a) is the union of four events associated to the four cases considered in Fig. [21 Each 
case in Eo involves the intersection with an outage event where the mutual information 
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between a user and the destination is below the rate R, except for case 4, where only 
the first frame is dedicated to user 1. 

b) follows directly from Q. 

c) uses the fact that the two frames in a block behave as parallel Gaussian channels whose 
capacities add together. Of course, both frames timeshare a time-interval, which gives 
a weight to each capacity term [10, Section 9.4], [39, Section 5.4.4]. 

([8]) follows from maximum ratio combining [39] at the destination during the second frame. 

■ 

The outage probability is obtained by integrating the joint probability distribution p(ai2, a2i, aid, a2d) 
over the volume defined by Eg. Just as for the Gaussian modulation, there is only one free 
parameter P because R and 7 are fixed by the protocol and the physical environment. Hence, 
given R and 7, one can optimize the value of p. For example, notice that for a low-SNR interuser 
channel, the outage probability improves while taking P smaller than 0.5 due to the enhanced 
protection of the source-relay channel. On the other hand, a P smaller than 0.5 results in lower 
achievable coding rates, as proved in Corollary [Tl The optimization of P, as already undertaken 
in [19] for Gaussian modulations, is not within the subject of this paper. 
There is an important conclusion to draw from the analysis of Prop. [TJ 

Corollary 1 In coded cooperation over a block-fading channel for the 2-user MAC with a 
cooperation level P, transmitting at a coding rate greater than min{P, {1 — P)) renders a single 
order diversity. 

Proof: It is sufficient to prove the stated corollary for coded cooperation over a Block 
Erasure Channel (BEC) [24], because a BEC is an extremal case of a block- fading channel. In a 
BEC, the fading gain a takes two possible values {0, +00}. An outage event on a point-to-point 
channel is defined by the fading gain a being zero. As a consequence, the possible values of 
the BPSK capacity on a BEC are confined to zero or one. Hence, for the two-user MAC, the 
mutual information Iid related to case 1 belongs to {l,P, (1 — /3),0}. A double diversity order 
is equivalent to stating that two outage events are necessary to lose the transmitted codeword. 
Take the scenario where the userl-to-destination channel has fading gain zero and the user2-to- 
destination channel has fading gain 00. In this scenario, the mutual information Iid is equal to 
p. All coding rates higher than P will limit the diversity order of the outage probability to one. 
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since only one channel in outage is enough to lose the codeword. From a similar reasoning, it is 
shown that Rc must be smaller than {1 — (3). This corollary is also valid for signaling strategies 
with M constellation points. ■ 

In the sequel, if not otherwise stated, we assume a rate equal io Rc= \- From Corollary [H 
we know that it is possible to have double diversity order for this rate if /5 G [|, |]. We stress 
on the fact that the proposed code construction is very flexible in parameters such as the block 
length and the coding rate, as long as this coding rate is smaller than min(/3, (1 — (3)). We will 
use (3 = 0.5 throughout this paper, which allows us to achieve diversity two for all rates Rc 
smaller than 1/2. We illustrate this in the numerical results by showing the WER performance 
of an LDPC code whose coding rate Rc approaches 1/2. 

IV. Full-diversity Low-Density coding for coded cooperation 

Codewords in coded cooperation are split over 2 frames. The first part of a codeword, 
transmitted during the first frame should protect information on the noisy source-relay channel. 
Consequently, a channel code, compatible with two distinct rates is to be devised. In non- 
cooperative communications, this property is known as rate-compatibility where parity bits of 
higher rate codes are embedded in those of lower rate codes [15]. This concept was heavily used 
in the context of adaptive coding and unequal error protection in packet data communication. 
The advantage is that all codes can be encoded/decoded using a single encoder/decoder which is 
for example an efficient framework for ARQ/FEC, where incremental parity bits are transmitted 
in response to erroneous decoding at the receiver. 

Rate-compatibility in the context of LDPC codes was first introduced by Li et al. [25] and Ha et 
al. [14] and further elaborated for example in [41]. Two techniques have been used: puncturing 
and extending. A fraction of parity bits of a mother code could be punctured to obtain higher 
rate codes. However, the resulting rate range is limited because the deletion of too many bits 
has a negative effect on decoding via belief propagation. To obtain a more dynamic range in 
rates, the technique of extending has been used. The extension is made by adding extra parity 
bits as illustrated in Fig. [3l where the overall code is the intersection of two constituent codes 
defined by H2 and Hi padded with zeros on the right. 

For simplicity, we only used the technique of extending to acquire rate-compatibility, but this 
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may be further optimized by combining puncturing and extending via known techniques [14], 
[25], [41]. 

A. Full-diversity LDPC codes 

Puncturing and erasing are slightly different in nature. From an algebraic decoding or a 
probabilistic decoding point of view, puncturing and erasing are identical, an erased/punctured bit 
is equivalent to an error with known location but unknown amplitude [26]. From a transmitter 
point of view, punctured bits have always fixed locations whereas channel erased bits have 
random locations. These similarities and differences in nature also exist for block erasures and 
block puncturing. 

Full-diversity coding on a relay channel must cope with block erasures. Consider the coding 
structure plotted in Fig. [3l If all parity bits 2p are erased due to deep fading in frame 2, then the 
decoder should be capable to retrieve information bits i thanks to Hi and possibly recompute 
2p thanks to if 2- Unfortunately, under deep fading in frame 1, a structure with a randomly 
generated H2 as in Fig. [3l cannot guarantee the retrieval of the information bits through H2. 
The aim of this section is to explain how H2 can be tuned in order to have full-diversity for 



Frame 1 



Frame 2 
2p 







Ip 



Hi 



H2 



Fig. 3. Parity-Check matrix of a rate-compatible LDPC code obtained by the extension of higher rate codes. Symbols are split 
into three classes: i for the information bits, Ip and 2p for two classes of parity bits. The classes i and Ip are transmitted by 
the source in frame 1. Parity bits 2p are transmitted in the second frame, for example by the relay after successful decoding of 
the first frame. 



March 9, 2009 



DRAFT 



13 



any left and right degree distribution and for any block length. The cases 2 and 4 considered in 
Fig. |2] always have double diversity order or more for user 1, because at least two outage events 
are needed to lose the codeword, one on the source-destination channel and one on the interuser 
channel. Since case 3 has always a better performance than case 1 for user 1, we will assume 
in the following the occurrence of case 1 where the transmission on the interuser channel in the 
first frame has been successful and both users are cooperating in the second frame. This means 
that, for the transmission of user 1, the first fading state of the overall codeword is given by the 
userl-to-destination channel and the second fading state of the overall codeword is given by the 
user2-to-destination channel. The key to full-diversity codes will be to spread the information 
bits of a user over both frames. 



Before describing the whole parity-check structure of a coded cooperation LDPC code in the 
next subsection, let us focus on the submatrix H2. We take the constituent code defined by H2 
to be a full-diversity LDPC code (referred to as root-LDPC code) as constructed by Boutros et 
al. in [9], [6] for non-cooperative single-antenna channels with two or more fading states per 
codeword. The Tanner graph notation for the root-LDPC code is given in Fig. SI This notation is 
essential for the analysis because we seek full-diversity under iterative decoding. Full-diversity 

' j Information node transmitted in frame 1, 
^ ^ connected to a rootcheck. 



Parity node transmitted in frame 1 . 



Parity node transmitted in frame 2. 



Information node transmitted in frame 2, 
connected to a rootcheck. 



Check node. 



Fig. 4. Notation for the Tanner graph of a full-diversity LDPC code. 






of a root-LDPC structure is created by rootchecks, a special type of checknodes in the Tanner 
graph. As shown in Fig. [5l the root and the leaves of this special checknode do not belong to the 
same frame. When the rootbit is in frame 1, the leavebits are in frame 2, and vice versa. Using 

March 9. 2009 DRAFT 



14 



the limiting case of a Block-Erasure Channel, it is easy to verify that a rootbit is determined 
via its rootcheck when its own frame is erased. The complete root-LDPC structure is built after 




Fig. 5. Two types of rootchecLs. On the left side, the rootbit belongs to frame 1 and the leavebits belong to frame 2. The 
symmetric case where channel states are switched is shown at the right side. 

splitting information bits into two classes, denoted li and 2i, and parity bits into two classes, 
denoted Ip and 2p. The checknodes are cut into two classes denoted 3c and 4^ The classes 3c 
and 4c consist of rootchecks for information bits li and 2i respectively. The complete root-LDPC 
structure including all types of nodes is illustrated in Figs. [6] and U\ Rootchecks are translated 
into two identity matrices (or permutation matrices in general) inside the parity-check matrix in 
Fig.lll 

The proof of full-diversity for block-Ray leigh fading can be found in [9]. Note that the diversity 
order of the root-LDPC code does not depend on the right or left degree distributions. For 
simplicity, we only showed a regular (3,6) structure in Fig. [6l 

For asymptotic code lengths, multiedge-type messages propagate in the root-LDPC graph [35]. 
One has to choose between two different root-LDPC ensembles. If we refer to the Tanner graph 
in Fig. [6l the two ensembles are distinguished as follows: (i) The first ensemble is built by 
two random edge permutations (edge interleavers) connecting 3c to (2i, 2p) and 4c to (li, Ip) 
respectively. This is equivalent to the random generation of two low-density matrices (H2i, H2p) 
and (Hii, Hip) in the parity-check matrix shown in Fig. U\ (ii) The second ensemble is built 
by four random edge permutations 3c — 2i, 3c — 2p, 4c — li, and 4c — Ip. In the root-LDPC 
parity-check matrix, this is equivalent to building seperately the four submatrices H2i, H2p, Hu, 
and Hip. For simplicity reasons, mainly in the density evolution (DE) analysis, we adopt the 

'The checknode notation Ic and 2c is reserved for Hi in the cooperative code as described in the next subsection. 
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Fig. 6. Tanner graph of a full-diversity LDPC code of length A'^ and rate i. This compact graph representation has been 
adopted from [9], [6], it is also known as protograph representation [38]. The integers labeling the edges of the Tanner graph 
indicate the degree of a node along those edges for a regular (3,6) root-LDPC code. The binary elements are split into four 
classes of each ^ bits. The checknodes are cut into two classes of ^ checks. 



Frame 1 Frame 2 



li 


ip 


2i 


2p 




1 

1 

1 









3c 


1 










Hii 




1 

1 

1 

1 





4c 



Fig. 7. Parity-Check matrix of a rate i root-LDPC code. 

first root-LDPC ensemble as part of the full-diversity cooperative code proposed in the next 
subsection. 

B. Rate-compatible full-diversity LDPC codes 

Now, the structure of a coded cooperation LDPC code is derived by joining the rate-compatibility 
property and the full-diversity property. In a first step, the global parity-check matrix is obtained 
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by embedding the root-LDPC matrix (Fig. |7]) into the rate-compatible matrix (Fig. |3]). This leads 
to an asymmetric code where class li may have have a higher coding gain than class 2i. In a 
second step, to get a balanced structure, we replace the zero-padded Hi by the direct sum of two 
rate Ri codes defined by His and Hir as illustrated in Fig. [8l Thus, the constituent code His 
protects bits li and Ip via extra parity bits p'l- Similarly, in the second frame, extra parity bits 
P2 are generated from 2i and 2p. The bottom of the global parity-check matrix simply includes 
the root-LDPC structure, connecting (li, Ip) to (2i, 2p). For simplicity we can assume that His 
and Hir belong to the same random rate Ri LDPC ensemble, defined by the degree distributions 
(Ai(x), Hence, if the degree distribution of the root-LDPC is {X2{x) , p2{x)) , we refer to 

the rate-compatible root-LDPC as a {\i{x), pi{x),\2{x), p2{x)) code. The Tanner graphs of a 
regular (3,9,3,6) LDPC code and an irregular (Ai(x), A2(x), p2(a;)) code are shown in 

Figs, m and [lOl Since we guarantee full-diversity via a root-LDPC with a fixed rate |, the global 
coding rate of the rate-compatible root-LDPC code observed at the destination is i?c = "y"- As 
a consequence, the global coding rate Rc can be easily varied through Ri. 







Frame 1 






Frame 2 






li 




ip 


p'l 


2i 


2p 










His 













Ic 











Hlr 


2c 


1 

1 

1 


1 








Hi, 


Hjp 





3c 




Hip 





1 

1 

1 

1 








4c 



Fig. 8. Parity-check matrix of a rate-compatible root-LDPC code for coded cooperation. The upper coding rate associated to 
His and Hir is -Ri — |, the bottom root-LDPC coding rate is i, and the overall coding rate is i?c = = |. 

Due to the identity matrices inside the parity-check matrix, new polynomials \2{x) appear in 
Fig. [To] in the connections li — Ac and 2i — 3c, as illustrated in Fig. [TTl 
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Proposition 2 In a Tanner graph with a left degree distribution \{x), isolating one edge per 
bitnode yields a new left degree distribution described by the polynomial X{x): 

i E,^i(j-i)/j 

Proof: Let us define Tbit,i as the number of edges connected to a bitnode of degree i. 
Similarly, the number of all edges is denoted Tbit. From section |Ill we know that A(x) = 
^4maa: \^rji.i-i exprcssBS the left degree distribution, where Aj is the fraction of all edges in the 
Tanner graph, connected to a bitnode of degree i. So finally Aj = -y^. A similar reasoning can 
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Fig. 10. Tanner graph of an irregular rate-compatible root-LDPC code for coded cooperation. The binary elements are split 
into six classes, p'l and p'2 of each and li, Ip, 2i, and 2p of each -^^j^ bits. The checknodes are cut into four 

classes of -^^^ checks. 




Fig. 11. Transition from a traditional representation, characterized by an edge distribution polynomial A (a;), towards a 
representation where one edge per bitnode is isolated resulting in a new degree distribution \{x). 



be followed to determine A, 



^, bit,i T ''i' 

^bit ~ X]j ^-^bit 

T T 
-ibit ~ Z^j ^^bit 



A, 

A,; 



V -i? — V 

j -J j 

H-1) 
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^) J2j 'f'^bit is equal to the number of edges that are removed which is equal to the 
number of bits. 

b) AjTbit is equal to the number of edges connected to a bit of degree i. 

■ 

In section |Vl we will also use p{x), which is defined similarly as A (x). 

Proposition 3 Consider a ( Ai (x) , pi (x) , A2 (x) , p2 (x) ) rate-compatible root-LDPC code for coded 
cooperation transmitted on a 2-user block-fading cooperative MAC. Then, under iterative belief 
propagation decoding, the rate-compatible root-LDPC code has full-diversity. 

Proof: As indicated in the design of full-diversity LDPC codes in subsection IIV-AI the 
diversity order of a root-LDPC code does not depend on its left or right degree distribution. 
Thus, we restrict this proof to a regular (3, 9, 3, 6) LDPC code. Let A^, i = 1 . . . — I denote 
the input log-ratio probabilistic messages to a checknode $ of degree dc- The output message 
A"^ for belief propagation is [33] 

A^=2th- (^nth(f)j, 

where th(x) denotes the hyperbolic-tangent function. Superscripts a and e stand for a priori and 
extrinsic, respectively. To simplify the proof, we show that the suboptimal min-sum decoder 
yields a diversity order 2. For a min-sum decoder, the output message produced by a checknode 
$ is now 

dc-l 

A^ = min(|A^|)nsign(An. 

i=l 

An information bit -d of class li has Aq = where Aq is the log-likelihood ratio coming 

from the likelihood p{ysr\'^). It also receives 6 messages A^^, z = 1 . . . 3 for its three neighbouring 
checknodes in each constituent code, j = 1,2 for the 2 constituent codes Hu and H2 respectively. 
The total a posteriori message corresponding to ^9 is A = Aq + X]j=i ZlLi Mj- The addition of 
Yl^=i-^ii cannot degrade the error probability Pe{li) because the convolution with the density 
of messages from His can only physically upgrade the resulting density. Thus, it is sufficient to 
prove that message Aq + J2t=i exhibits full-diversity, which is proven in [9]. ■ 
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V. Density Evolution on the Block-Fading Relay Channel 

Richardson and Urbanke [32], [33] established that, if the block length is large enough, 
(almost) all codes in an ensemble of codej^ behave alike, so the determination of the average 
behavior is sufficient to characterize a particular code behavior. This average behavior converges 
to the cycle-free case if the block length augments and it can be found in a deterministic way 
through density evolution (DE). The evolution trees represent the local neighborhood of a bitnode 
in an infinite length code whose graph has no cycles, hence incoming messages to every node 
are independent. 

A. Interuser channel 

To determine the density of messages propagating in the graph of the constituent code His, 
the following notation is used: 

'^^(^) = density of message from a bitnode to a checknode in the m* iteration. (11) 
^sr{x) = density of the likelihood of the source-relay channel. (12) 

Let Xi ~ pi{x) and X2 ~ pi{x) be two independent real random variables. The density function 
of Xi + X2 is obtained by convolving the two original densities, written as pi(x) ®p2{x). The 
notation p(a;)'^" denotes the convolution of p{x) with itself n times. 

Let Xi ~ pi{x) and X2 ~ P2{x) be two independent real random variables. The density function 
p{y) of the variable Y = 2 th~^ (th th (^)), obtained through a checknode with Xi and X2 
at the input, is obtained through the R-convolution [33], written as pi{x) Qp2{x). The notation 
p(x)®'^ denotes the R-convolution of p{x) with itself n times. 

To simplify the notations, we use the following definitions: 

''The ensemble of all LDPC -codes that satisfy the left degree distribution A(a:;) and right degree distribution p{x) is considered. 
The ensemble is equipped with a uniform probability distribution. 
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A {p (x)) = 5^ A. pixr'-\ pipix)) = J2 P^ Pixf'-'. 

i i 

In the next subsection we will also use the following definitions: 

p{p{x),t{x)) = Y,{p,p{xf'-^Qt{x)), 

i 

\*{p{x)) = A(p(x))®(p(x)), P* {p (x)) = p {p (x)) Q {p (x)) . 

The first definition is necessary because of the non-linearity of the R-convolution. Therefore, the 
first equation is not equal to t(x) & p (p (x)). The next subsection will also use the polynomials 
p* (x) and A* (x) which are defined by combining the two transformations, denoted by (.) (see 
introduction) and (.)*. 

Fig. [T2l illustrates the local neighborhood of a bitnode in the constituent code His. 



Ic 




Fig. 12. Local neighborhood of a bitnode in the constituent code His- This tree is used to determine the evolution of density 
dsrix) of messages from a bitnode to a checknode. 



The DE equation in the neighborhood of the bitnode for a (Ai(x), LDPC code [32] is. 
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for all m, 

=/i,,(x)® Ai(^Pi(C(a;)))- (13) 

The threshold of a code is the minimum SNR at which a codeword can be decoded perfectly 
[32]. Comparing the received signal-to-noise ratio with this threshold, the relay and the source 
can determine whether the interuser transmissions can be decoded successfully and consequently 
decide what to transmit in the second frame. 



B. Overall two-user MAC 

The proposed {Xi{x) , pi{x) , X2{x) , p2{x)) root-LDPC code has 6 variable node types and 4 
checknode types. Consequently, the evolution of message densities under iterative decoding has 
to be described through multiple evolution trees for a binary root-LDPC code. Figs. [131 [El [HI 
[TtI [TSl and \T9\ show the local neighborhood of classes li, Ip, and p[. The local neighborhood 
of classes 2i, 2p, and P2 are equivalent because of code symmetry. 

To determine the density of messages, the following notation is used: 









[x) 


= density of message from li to Ic and 2i to 2c respectively, at the m"^ 


iteration, 


/r( 


x). 




[x) 


= density of message from li to 3c and 2i to 4c respectively, at the m* 


iteration, 








[x) 


= density of message from li to 4c and 2i to 3c respectively, at the m"^ 


iteration, 




»' 




[x) 


= density of message from Ip to Ic and 2p to 2c respectively, at the m* 


' iteration 


/r 


(x) 




{x) 


= density of message from Ip to 4c and 2p to 3c respectively, at the m'*^ 


' iteration, 




[x).. 


,92"! 


{x) 


= density of message from p[ to Ic and to 2c respectively in the m'*^ 


iteration, 








[x] 


= density of the likelihood of the channel in the i'th frame. 





Note that fi2{x) depends on the success or the failure of the transmissions in the first frame. 
In Figs. \T3\ \T5\ [Ml the local neighborhood of a bitnode of the class li is given. 

Proposition 4 The DE equations in the neighborhood of li, Ip and p'l for a ( Ai (x) , pi (x) , A2 (x) , p2 (x) ) 
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Fig. 13. Local neighborhood of bitnode li. This tree is used to determine the evolution of the density of messages li — > Ic. 



root-LDPC code for coded cooperation are, for all m, 



aT^\x) = fiiix)^X2\^P2{fu,c9Ti^) + fip^clTi^),fTi^)) 
®Ai (^pi (/Hie aT{x) + /ipie kT{x) + f,,,, qT{x)) 

®f>2 (/2.3c 9"2\^) + /2p3c 1^{X)) , (14) 
®A2fp2(/H4c^7r(^) + /lp4c/r(^), f^i^))]. (15) 
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®A2(p2(/i.4c gTi^) + hpAc fT{x)) 

®P2 (fe 92 {X) + /2p3c I'^ix)) , (16) 

= Pi{x)®X;(^p2{fu4c9Ti^) + fmclTi^),fTi^)) 

^\i(pi{hic «r(x) + + gr(a;))J> (17) 

lT^\x) = ^^l{x)®\l(^p^{f^acaT{x) + ^p,,kT{x) + fp,^,qT{x)) 

®X2(^P~2{fu,cgTi^) + fip^clTi^), fTi^))^, (18) 
qT^\x) = p.iix)0\i(pi{hicaT{x) + fipickT{x) + fp>i,qTix))) (19) 



where 

/i.4. = (20) 

/l.4c = (21) 

/lilc = /iplc, (23) 

fp[lc = 1 ~ /lilc — /iplc, (24) 

/2i3c = /li4c, (25) 

/2p3c = /lp4c- (26) 

Proof: Equations (fT4l)-(|26l) are directly derived from the local neighborhood trees. To obtain 
the proportionality factors (|20l)-(|26l). it is important to remark that we use the first ensemble of 
root-LDPC codes, as explained at the end of section IIV-A[ Let T denote the total number of 
edges between the variable nodes (li — lp) and the checknodes 4c. Fig. [14] illustrates how fip^^ 
and /ij4c are obtained: 



March 9, 2009 



DRAFT 



25 




4c 



Fig. 14. Part of the compact graph representation of the Tanner graph of a root-LDPC for coded cooperation. The number of 
edges connecting (li, Ip) to 4c is T. the number of edges connecting Ip to 4c is Tip. The number of edges connecting li to 
4c is Tu. 



T 



a) RiN/A 



Tip 

flpAc 
fliic 



a) RiN/A 
a) RiN/A 

VI Tip 
T 

6) Tu 
T ■ 



(27) 

(28) 

(29) 

(30) 
(31) 



a) The number of checknodes connected to i edges of T is ^T. A Similar reasoning 
proves equations (l28l) and (|29l ). 

b) The fraction of edges T connecting Ip to 4c is fip^c- The fraction of edges T connecting 
li to 4c is /ii4c. 
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Fig. 16. Local neighborhood of bitnode li. This tree is used to determine the evolution of the density of messages li 4c. 
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Fig. 18. Local neighborhood of bitnode Ip. This tree is used to determine the evolution of the density of messages Ip — > 4c. 
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VI. Numerical Results 

A. Infinite Length Codes 

We evaluated the asymptotic performance of rate-compatible root-LDPC codes by applying 
DE on the proposed code construction. We studied different scenarios: 

1) Scenario 1: 

• The average SNR of the independent interuser channels is 5dB higher than the average 
SNR on the source-destination link. 

• The average SNR of the relay-destination link is equal to that on the source-destination 
link. 

• The coding rate is _Rc = | and the cooperation level is /3 = 0.5. 

Fig. [20] shows the main results: the word error rate (WER) of a regular (3,9,3,6) rate-compatible 
root-LDPC code and of an irregular (Ai(x), pi(x), X2{x) ■, P2{x)) rate-compatible root-LDPC code 
with left and right degree distributions given by the polynomials 

Ai(x) = 0.1989X + 0.2305x2 + 0.0068x^ + Q.277Ax^ + 0.14267a;^^ + 0.1335x2° + 0.0102x2\ 

Pi{x) = x^^ 

A2(x) = 0.22767X + 0.20333x2 + 0.2145x^ + 0.011048x^ + 0.34346x^^ 
P2(x) = 0.5x^ + 0.5x1 

As can be seen clearly on the figure, the regular rate-compatible root-LDPC code is always 
within L5dB from the outage probability limit, whereas the irregular rate-compatible root-LDPC 
code is always within IdB from the outage probability limit. This distance is respected for 
many variations of the channel conditions, such as other interuser channel conditions or uplink 
channel conditions. Note that our code construction can, just as in [16], be applied on a fuU- 
uplex channel, doubling the overall spectral efficiency. As mentioned before, the coding rate is 
adjustable by varying the number of parity bits p[ and p'2, which is illustrated in scenario 2. 
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2) Scenario 2: 

• The average SNR of the independent interuser channels is 12dB higher than the average 
SNR on the source-destination link. 

• The average SNR of the relay-destination link is 4dB higher than the average SNR on the 
source-destination link. 

• The coding rate is Rc = 0.45 and the cooperation level is /3 = 0.5. 

Here, we imitated the channel conditions used in [16]. These channel conditions are optimistic, 
allowing a high coding-rate for the source-relay channel. Fig. [21] illustrates the WER of an 
irregular (Ai(x), pi(x), X2{x), p2{x)) rate-compatible root-LDPC code with left and right degree 
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distributions given by the polynomials 



Ai(x) 

p2{x) 



0.1581a; + 0.2648a;^ + 0.1116a;^ + 0.1354a;'^ + 0.3301a;^^ 



X 



43 



0.234413X + 0.21392x2 + 0.123711a;^ + 0.125548x^ + 0.30241x^^ 
0.71875x^ + 0.28125x1 



The coding rate for the interuser channel subcode Hi is equal to 0.9. 



LU 
T3 
O 




12 15 
E3/N0 (dB) 

Fig. 21. Density evolution of rate-compatible root-LDPC codes with coding rate Rc — 0.45 with iterative decoding on a 
cooperative MAC with two users. Es/Nq is the average symbol energy-to-noise ratio on the source-destation link. 



Fig. |2l] shows that a rate-compatible root-LDPC code with a coding rate approaching | still 
exhibits a near-outage performance. A code is referred to as MDS code [20], [27] when the 
diversity order is equal to [Nu{l — Rc)\ + 1, where Nu is the number of users. In contrast to 
code constructions previously published in the literature where no root-LDPC codes are used, 
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this code construction renders an MDS code (the diversity order is still two if a coding rate 
Rc = 0.5 is adopted), because the diversity order only depends on the constituent code H2, 
which has been shown to be MDS in [9]. However, it is not possible to show DE results for the 
rate-compatible root-LDPC with a coding rate Rc = 0.5, because this implies that Ri = 1 and 
that there is no threshold effect on the interuser channel. 

B. Finite Length LDPC Codes 

It is interesting to evaluate the finite length performance of the proposed rate-compatibel root- 
LDPC codes. Not only to approve the asymptotic performance, but also to see how to generate 
an instance of the parity-check matrix, given by Fig. [8l Before showing the results, we will first 
discuss the practical generation of this parity-check matrix. 

Consider case 1 from Fig. [2l For the decoding process, the destination will apply the sum-product 
algorithm on the overall graph including His, Hir, and H2. For the encoding process, it is easier 
to determine the parity bits p[, p'2, and {lp,2p) with the parity-check matrices His, Hir, and 
H2 respectively. As with standard LDPC encoding, these matrices will then be systemized to 
determine the parity bits. An important constraint for the decoding process is the alignment in the 
overall parity-check matrix of common bit nodes in both constituent codes. This can be achieved 
by prohibiting column permutations during the systemization of His, Hir and H2. Except for 
case 4, which only decodes on H^, the other cases need the same constraints. 

1) Generation of Hu and Hir: His and Hir are randomly generated satisfying the degree 
distribution pi{x) for its rows and the degree distribution Ai(x) for its columns. A sufficient 
condition to prohibit column permutations during the systemization of Hu and Hir is imposing 
on Hpi^ and Hpi^ to be full-rank. Hpi^ (Hp'^ respectively) is the most right square matrix of His 
(Hir respectively). 

2) Generation of H2: The generation of H2 can be split in the generation of H^c and H^c, 
where H-^c {H^c resp.) is the upper part (resp. lower part) of the parity-check matrix H2- H^c 
is the concatenation of an identity matrix (permutation matrix), zeros and a randomly generated 
matrix {H2i,H2p). The rows of {H2i, H2p) satisfy the degree distribution p2{x), the columns of 
the most left square matrix H2i satisfy the degree distribution A2(x) and the columns of the most 
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right square matrix satisfy the degree distribution A2(x). This is equivalent to generating a 
random graph with two classes of bitnodes at the left side and one class of checknodes at the 
right side of the graph. If 77,3c is the number of checknodes at the right side, then a random graph 
with J^'-^- edges is generated. A fraction of the edges is connected to bit nodes of the 

class 2i, whereas a fraction of the edges is connected to bit nodes of the class 2p. In 

the end, the identity matrix is simply added. H^^c is generated similarly. 

For the encoding process, we have to systemize this matrix. One solution is to switch the 
columns associated with the li bit node class and the 2p bit node class. The most left square 
matrix of H2 will then be block-diagonal with and Hip on its diagonal. Having H2p and 
Hip full-rank is consequently a sufficient condition to exclude column permutations during the 
systemization of this matrix. After the generation of (2p, Ip), all the bits are put in the required 
order li — Ip — 2i — 2p by switching back the bits of the classes li and Ip. 

3) WER performance of finite length LDPC codes: The WER performance of different block 
lengths is illustrated by Fig. |22l We used the same channel conditions and coding rate as in 
scenario 1 in the previous subsection. This figure approves the DE results for an asymptotic 
block length. 

C. Comparison with previous work 

As mentioned in the introduction, especially rate-compatible punctured convolutional codes 
(RCPC) have been used in coded cooperation. The main drawback of these codes is that the WER 
increases with the logarithm of the block length to the power d where d is the diversity order 
[7], [8], whereas the WER of near-outage codes should be independent of the block length. This 
can be seen clearly on Fig. [23l where we show the WER of two rate-compatible non-recursive 
non-systematic (75,53,47) convolutional codes with blocklength 500 and 5000 respectively. We 
used the same channel conditions and coding rate as in scenario 1 in subsection IVI-AI 
We also compared with another protocol. Decode and Forward (DF), using near-outage LDPC 
codes for this protocol. Despite the fact that this implementation has near-outage performance, 
the WER performance is worse than that of our code construction. The reason is that the outage 
probability limit of DF is higher than that of coded cooperation. 
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Fig. 22. Comparison of rate-compatible root-LDPC codes for different block lengths with iterative decoding on a cooperative 
MAC for two users, coding rate Rc — 1/3. The ratio Et/No is the average information bit energy-to-noise ratio on the 
source-destation link. 
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Fig. 23. Comparison of rate-compatible root-LDPC codes for coded cooperation with other work on a cooperative MAC for two 
users. We simulated LDPC codes for Decode and Forward under iterative decoding and an implementation of rate-compatible 
convolutional codes [17]. The ratio Eb/No is the average information bit energy-to-noise ratio on the source-destation link. 



March 9, 2009 



DRAFT 



36 



VII. Conclusion 

We have studied LDPC codes for relay channels in a slowly varying fading environment 
under iterative decoding. We have introduced the new family of rate-compatible root-LDPC 
codes, which combines the rate-compatibility property with the full-diversity property for any 
coding rate Rc < Rcmax = min(/5, 1 — (3), where (3 is the cooperation level. We have shown that 
the error rate performance of regular and irregular rate-compatible root-LDPC codes is close to 
the outage probability limit and this occurs for all block lengths (finite and infinite) and all rates 
not exceeding Rcmax- Its flexibility and high performance makes rate-compatible root-LDPC 
attractive for wireless cooperative communications scenarios with slowly varying fading. 
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